Three-dimensional static vortex solitons in incommensurate magnetic crystals 
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A new type of three-dimensional magnetic soliton in easy-axis ferromagnets is predicted by taking 
simultaneous account of the Dzyaloshinsky-Moriya interaction and an external magnetic held. The 
numerically obtained static three-dimensional solitons with a finite energy are characterized by a 
Hopf topological index H = and have a vortex structure. The structure of these solitons and the 
dependence of their energy on the external held are determined. The asymptotic behavior of these 
solitons is investigated and a necessary condition for their existence is found. 



I. INTRODUCTION 

At present a large number of magnetically ordered 
crystals without an inversion center are known in which 
exchange-relativistic interactions lead to the formation 
of long-period magnetic structures whose periods are in- 
commensurate with the crystal-chemical periods |l], ^). 
In the energy for this system the exchange relativistic 
interaction is described by terms which are linear in the 
first derivatives of the magnetization (the Lifshitz invari- 
ant). The possibility that they might have a key role in 
the formation of helicoidal magnetic structures was first 
advanced in || . Analytic solutions || showed that for 
certain values of the Dzyaloshinsky constant modulated 
structures can exist in the system with incommensurate 
periods. It was shown later || |J that twodimensional 
vortices can exist against the background of an incom- 
mensurate structure in magnets with an "easy-plane" 
anisotropy according to a two-dimensional generalization 
of the Dzyaloshinsky model (2D sine-Gordon). It has 
been shown j^, ^) that two-dimensional magnetic vortices 
can exist in magnetically ordered crystals with an "easy 
axis" anisotropy for a certain range of external magnetic 
fields. Depending on the parameters of the ferromagnetic 
medium and the external magnetic field, both unstable 
and stable vortex states with finite energies can develop 
fio| , |i~t| . The structure of threedimensional solitons and 
their domains of existence in these magnetic materials 
have not yet been studied. 

It is known that the Landau-Lifshitz equation without 
dissipation, 



dM 



- 7 [M X H eff ] , H eff 



SE 
"SM' 



7 >0 (1) 



allows solutions in the form of three-dimensional proces- 
sional solitons with a stationary profile under the fol- 
lowing conditions. If the energy E of the ferromagnet 
includes only the exchange energy, 



E exch = (a/2) / (dMY dr 



(2) 



i.e., the ferromagnet is isotropic, then |12| nonstationary 
solitons can exist which move uniformly along the axis 



of anisotropy at a nonzero velocity. Here, however, the 
energy is proportional to the linear size of the soliton and 
the question of stability with respect to collapse is espe- 
cially acute. In addition, it is not entirely right to neglect 
the energy of magnetic dipole interactions, since demag- 
netizing fields can facilitate the collapse of a soliton. 

If, on the other hand, the energy E of the ferromagnet 
includes the uniaxial magnetic anisotropy energy, 



= 08/2) / (M x 2 + M y 2 ) dr, p > 0, (3) 
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as well as the exchange energy E exc h then solitons with 
a different structure can develop, and which either move 
uniformly jl3[ 14 or are stationary |l5|-|l8|. The lin- 
ear sizes of these solitons are no longer determined by 
the value of an arbitrary parameter, but depend on the 
characteristic magnetic length Iq = y/a/j3. The contri- 
bution of the magnetic dipole interaction energy is really 
insignificant here when the magnetic medium has a high 
quality factor Q = /?/ (4tt) . 

Nevertheless, in all the above cases we are speaking of 
dynamic solitons, i.e., of the case where, when both the 
precession frequency uj = and the velocity V = 0, a 
soliton cannot exist. This has seriously inhibited exper- 
imental studies. A real ferromagnetic medium is always 
dissipative and the lifetime of a dynamic soliton may be 
very short. In order to slow down or entirely stop the 
dissipation of solitons, it is necessary to apply a special 
external actuation to the medium. In addition, the nec- 
essary condition for observing solitons continues to be 
the development of some mechanism for their rapid for- 
mation, so that the average lifetime of a soliton exceeds 
its formation time |]l6| . 

Static three-dimensional solitons cannot exist in any of 
the above models of ferromagnets. This follows immedi- 
ately from the Hobart-Derrick theorem 19, 23], which is 
also known as the virial theorem plj . 

In many cubic ferromagnets without an inversion cen- 
ter, the Dzyaloshinsky-Moriya energy makes a significant 
contribution to their energy 122H24 



Edm = 



D J M ■ (V x M) dr, 
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where D is the Dzyaloshinsky constant. Then the total 
energy of the ferromagnet is 



Er 



(5) 



axially symmetric distribution of the polar angle of the 
magnetization and a vortex structure for the azimuthal 
angle 



where Eh is the energy in an external magnetic field of 
strength Hp 6 -* = (0, 0, H ) relative to the ground magne- 
tization state M. in f = (0, 0, M ): 



Q = 9(r,z), $ = if + cj)(r,z) 
with the boundary conditions 

9^0 (|r|->oo) 



Eh = Hq J (Mo — M z ) dr. (6) at infinity and the symmetries 



It turns out p5fl that in the three-dimensional case 
the functional (19) does not fall under the prohibition of 
the Hobart-Dcrrick theorem, so the existence of local- 
ized states remains an open question. Here the necessary 
condition is Edm < p6| . 

This paper is a study of static three-dimensional soli- 
tons in the model of (pj). It is organized as follows: The 
structure of the solitons, their stability, and the magnetic 
field dependence of their energy are discussed in Section 
II. In Section III the asymptotic behavior of the solutions 
is studied. The resulting solitons are analyzed topologi- 
cally in Section IV. 



II. THREE-DIMENSIONAL STATIC VORTEX 
SOLITONS 

It is reasonable to assume that three-dimensional soli- 
tons, like plane vortices, can be in stable or unstable 
states. The latter cannot be found by direct methods of 
unconditional minimization of the energy functional. In 
order to avoid this restriction, we shall seek a solution 
of the auxiliary variational problem of minimizing the 
functional 



Edm, 



with a fixed value of the integral 



N = J (Mo - M z )dv. 



(7) 



(8) 



We parametrize the magnetization vector M in terms of 
the angular variables 9 and $, as 



M = Mo ■ m = Mo ■ (sin®cos$, sinQsin®, cosQ) (9) 

In the presence of (|J) , the energy density of the ferro- 
magnet is invariant with respect to simultaneous spatial 
rotations by an arbitrary angle e and the spatial coordi- 
nates, and the magnetization vector, to 



ip — > ip + e, $ ^ $ 



(10) 



where if is the polar angle of a cylindrical coordinate 
system (r, ip, z). Thus, we shall study solitons with an 



0(r, 



(r,z), 



(r,z). 



(11) 



(12) 



(13) 



No boundary conditions arc imposed initially at the axis 
of symmetry but show up as a result of the numerical 
procedure. 

The functions $ and 9 are solutions of the posed varia- 
tion problem and, at the same time, extremal functionals 
of the energy (pj) for a certain value of Ho- The problem 
actually involves finding the functions (f>(r, z) and 9(r, z) 
that specify the unit vector field on the half plane 



n(r, z) = (sinOcosrf), sin8sin<fi, cos9), 



(14) 



To calculate the value of Ho corresponding to extremum 
of the energy functional, we use the necessary condition 
for an extremum: 



—E(m(r + A • r)) 
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from which we immediately obtain 
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The variational problem was solved by the same nu- 
merical method as for conditional minimization in [ [l8[ . 
The initial configuration for the field was specified by 
the same functions as in 113 ] for determining the struc- 
ture of moving magnon droplets. To test the resulting 
configurations, we calculated the corresponding external 
magnetic field strength Ho according to (|l6|), and, at the 
grid points, the values of the parametrizing angles (j>(r, z) 
and 9(r, z) and their first and second order spatial deriva- 
tives. After this, the discrepancy in the Euler-Lagrange 
equations for the energy functional (0) was calculated 
directly: 
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FIG. 1: A x~ h diagram of the results of the numerical 

calculations. Each point corresponds to a successful 
computational cycle, i.e., to establishing the existence of 
a soliton and determining its structure. 



where the differential operator 

dr 2 + dz 2 
and the dimensionless parameters 
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as in jLl|. 

In calculations using the method described here it is 
necessary every time to specify the value of the integral 
(pj) and the parameter \. And only if the calculation 
leads to a positive result (to a soliton solution) can be 
calculate a corresponding value of the parameter h. In 
addition, if minimizing the functional does not yield a 
soliton solution, this does not mean the latter is nonexis- 
tent. This sort of situation can arise, for example, when 
the modelling region is not the right size (the soliton is 
too big). In light of this, only points corresponding to 
positive results are plotted in Fig.[l]. 

Figure || shows contours of constant angle 9 for a typ- 
ical soliton in the (r, z) half plane for two values of the 
external magnetic field. In a three-dimensional space a 
constant value of the angle 6 corresponds to one or sev- 
eral toroids with a complicated envelope. 

The profile 9(r, z) of a soliton for z = shown in Fig.|| 
shows that the characteristic size of the localized struc- 
tures being studied here is on the order of 10Zo- 

Figure [| is a plot of the energy of the soliton as a 
function of the external magnetic field. As the magnetic 
field is increased, the energy also rises. As these graphs 
show, to excite a soliton with a given energy in a material 
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FIG. 2: Contours of the angle 9 of a soliton 
configuration for two values of the parameter h with 
constant X = 1-96: h = 0.60 (a), h = 0.62 (b). 



with a higher value of D, it is necessary to increase the 
external magnetic field strength Hq. 

It turned out that all the solitons that were found are 
stable with respect to scaling perturbations: 



— £(m(r+Ar)) 



> 0, 



A=0 



(21) 
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FIG. 3: 9 as a function of the coordinate r at z = for 
a soliton with the parameters of a ferromagnetic 
medium x = 1-57, h = 0.25. 



200 



3 



150 



200 



150 



100 



100 1 1 1 

0.22 0.24 0.26 0.28 0.58 0.6 0.62 0.64 

h h 

(a) (b) 

FIG. 4: Soliton energy as a function of the parameter h 
for two values of the parameter \ for a ferromagnetic 
medium: 1.571 (a) and 1.963 (b). The points indicate 
the results of numerical calculations. 
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Nevertheless, the extrema that were found do not provide 
a minimum for the energy functional (^|) in a strict math- 
ematical sense. We examined the stability of the solitons 
with respect to small perturbations. A special type of 



functions 5<fi(r, z) and S9(r, z) (complex perturbations) 
was found such that 



7 E(<f) + \6<f>,9+ X69) 



<0. 



A=0 



(25) 



However, even if the structures being studied here turn 
out to be unstable, they may have a rather long lifetime. 



III. ASYMPTOTIC BEHAVIOR 

Data from the numerical calculations show that with 
distance from the center of a soliton the azimuthal angle 
of the vector n (|l4|) becomes a linear function of the 
coordinate z, i.e. (j)(r, z) —¥ lixm + 7r/2 + kz as |r| — > oo, 
m G Z. In addition, the soliton structure found as a 
result of minimization shows that 9(0, z) — 0. Since the 
angle 9 goes monotonically to zero as |r| — > oo, the initial 
system of Eqs. ( |l7| ) and ([l8]) can be linearized for small 
9 with larger values of |r|: 



4xk 
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9 + A 2 9 = 



(26) 



(-(4x/7rf o )+2fc)0,0 = O. 



(27) 

It follows at once from (|2^) that the constant k is deter- 
mined by the Dzyaloshinsky constant, as 

k = 2x/(vr/ ) = D/a. (28) 

Using the auxiliary function 

tp(r, z, ip) = 9(r, z) cos ((f) (29) 

it is convenient to rewrite Eq. ( |2^ ) in the form 

= 7 2 V>, (30) 

where A is the laplacian operator in three dimensional 
space and 



7 = \/l + /l-(4 X 2 /7T 2 )A). 



(31) 



We shall seek a solution for (|29|) and ( |30| ) subject to 
the conditions 



f^^0,r^0, 

9(0, z) = 0, 

9^0 (|r|^oo), 
[9(r,-z) = 9(r,z). 



(32) 



Note that in the special case of 7 = 0, Eq. (pQ) trans- 
forms to the Laplace equation 



A^ = 0, 



(33) 
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and has an interesting electrostatic analog: the function 
if> is the potential for the electric field of a system con- 
sisting of charges distributed along a circle of arbitrary, 
but finite, radius r — a, at z = 0, with a linear density 
proportional to cos(ip) and an infinitely long, grounded 
filament passing along the axis of symmetry, where the 
potential of the filament, like the potential at infinity, is 
zero. In fact, Eq.(p3|) permits separation of variables in 
toroidal coordinates (it, v,ip), 

tanh(u) = (2ar)/(a 2 + r 2 + z 2 ), < u < oo, (34) 
cot{v) = (a 2 - r 2 - z 2 )/(2az), < v < 2ir, (35) 

and has the exact solution lETf : 



i/j = C^/a/rQ 1/2 (x)cos(ip), 



(36) 



where C is an arbitrary positive constant, \ = coth(u), 
Qi/2(x) is the Legendre function of the second kind of 
order 1/2. This function can be written as the integral 



Qi /2 (x) 
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l+X 
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where K and E are the complete elliptical integrals of the 
first and second kind, respectively. The solution (|36|), like 
every linear combination of similar solutions, satisfies Eq. 
( |32"| ) and yields the following asymptotic representation 
for the function 



(r,z) ~ (r 2 + z 2 ) 



-3/2, 



(M 



3), 7 = 0. (38) 



The result (38) can be derived formally in another way. 



If d is the polar angle of the spherical coordinate sys- 

- z 2 , then separating 
yields a basis of the 



tern (p,-&,tp) with p = |r| = vr 
the variables in (|3^) and using 
functions of the form 



i' n = Cp-^P^COS^COsiifi). 



(39) 

Given the imposed conditions ([32]) , n is a positive odd 
number and in the asymptotic limit, only n = 1 "sur- 
vives" . This immediately yields the asymptote ( 38 ) . 

If 7 7^ 0, then the Helmholtz equation ( |30| ) is solved by 
separation of variables in spherical coordinates. Then, 
given (p9|), we find the basis of the functions in the form 

C 

ip n = — K n+1/2 (jp)Pn(cosD)cos(<p), (40) 

where K v (x) is the modified Bessel function of the second 
kind. As a result, we obtain the final formula for the 
asymptotic behavior: 



9(r,z) 



(1 + 7V r 2 + z 2 ) r 

e ~/Vr' 2 +z 2 ( r 2 _|_ z 2)3/2 



(|r|^oo). (41) 



This implies the necessary condition for the existence of a 
three-dimensional vortex soliton: 7 > 0. Then, according 
to ©, 



1 + h - (4 X 2 /tt 2 ) > 0. 



(42) 



This result is in qualitative agreement with the plot in 
Fig.|l|. Note that this condition is the instability criterion 
for a helicoidal structure in uniaxial magnets |J . 
When the condition 



1 + h - (4x 2 /tt 2 ) < 0, 



(43) 



holds, the ground state of the hamiltonian (||) corre- 
sponds to the helicoidal structure 



cos(e) 



(4 X 2 /tt 2 )-1 



$ = kz. 



(44) 



where the step size of the spiral A = 27r/fc and k is the 
same as in (p8f). This sort of structure can be classified 
as a ferromagnetic spiral (FS) [Q, |j. Here the energy 
density 



AV/3(l + /i-(4 X 2 /7r 2 )) 2 
2tt 2 (4x 2 /tt 2 )-1 



w 



(45) 



This implies that the condition (^2|) holds, a uniform sta- 
ble state 8 — is realized in the system. Thus, the soli- 
tons we have been studying are localized excitations of a 
uniform state. 

The volume energy density of the soliton, vo, falls off 
exponentially with distance from the center of the soliton, 



w(p,tf) 



cos(kp cos"d){sindY 
e~<Pp 



(p^oo). (46) 



IV. THE TOPOLOGY OF VORTEX SOLITONS 

The unit vector field m is continuous and specified at 
each point of K 3 space, while in the asymptotic limit 
m — > (0, 0, 1) as |r| — > 00. These kinds of fields map the 
space K 3 U {00} onto a two-dimensional sphere S 2 are 
classified as belonging to homotopic class ^(S 2 ) = Z, 
and are characterized by an integral Hopf index H: 



H = 



1 



(87T 



F • Arfr, 



(47) 



where Fi — e^fcrn • (V,-m x Vfcm) and (V x A) = 2F. 
The expression for the Hopf index H of the fields ( |ll| ) 
simplifies pa, B(j to 



H=±- 

4ir 



00 Jo 



n • \d r n x d 7 vi\drdz 



(48) 



If the index H in (|48|) is nonzero, then the soliton can 

-|33 



be classified as a toroidal hopfion |14j, |29|-|33j . A numer- 
ical calculation shows that for all the solitons obtained 
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FIG. 5: The distribution of the magnetization for a 
soliton with h — 0.11, x — 1-37. The toroids 
constructed in the interval < tp < tt correspond to two 
values of the angle 0: for the outer toroid 9 = 0.3 and 
for the inner, 8 = 2.0 



here, H = 0. An example of a nontopological soliton 
with H — is a vortexfree magnon droplet . For the 
class of objects studied here, as in Fig.|[ one can see a 
structure of vortex wheels. The solitons found here can 
be classified as topological in the sense that they form an 
ensemble of "hopfion-antihopfion" pairs. 



V. CONCLUSION 

It has been shown here that, within a definite region of 
parameters, three-dimensional vortex solitons can exist 
in incommensurate magnetically ordered crystals. The 
structure of topological solitons with a nonzero Hopf in- 



dex H (hopfions) will be described in later papers. The 
static solitons obtained here do not have the disadvan- 
tages of precessing magnon droplets [pi - Thus, despite 
their possible instability, they may have a sufficiently 
long lifetime that they can be observed in real exper- 
iments. In particular, if the external magnetic field is 
changed suddenly to a state ( [42] ) for a stable helicoidal 
structure (^4j), then the magnet will end up in an un- 
stable phase and the transition to a uniform state may 
be accompanied by the appearance of solitons. It is also 
possible that within some range of the parameters x ~ h 
that has not been studied in this paper, solitons of this 
type will be stable from a mathematical standpoint. 

We thank B. A. Ivanov for the invitation to submit this 
article for the issue of this journal devoted to the 80-th 
birthday of V. G. Baryakhtar, whose scientific activity 
and school have always made significant contributions to 
the development of the physics of magnetic phenomena. 
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